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The ferrimagnetic spinel oxide Zn^Fea-^CU combines high Curie temperature and spin polariza- 
tion with tunable electrical and magnetic properties, making it a promising functional material for 
spintronic devices. We have grown epitaxial Zn^Fea-zCU thin films (0 < x < 0.9) on MgO(OOl) 
substrates with excellent structural properties both in pure Ar atmosphere and an Ar/C>2 mixture 
by laser molecular beam epitaxy and systematically studied their structural, magnetotransport and 
magnetic properties. We find that the electrical conductivity and the saturation magnetization can 
be tuned over a wide range (f0 2 . . . f0 4 _1 m _1 and f .0 . . . 3.2 /xe/f-U. at room temperature) by Zn 
substitution and/or finite oxygen partial pressure during growth. Our extensive characterization of 
the films provides a clear picture of the underlying physics of the spinel ferrimagnet Zn I Fe3_ a: 04 
with antiparallel Fe moments on the A and B sublattice: (i) Zn substitution removes both Fe^ + 
moments from the A sublattice and itinerant charge carriers from the B sublattice, (ii) growth in 
finite oxygen partial pressure generates Fe vacancies on the B sublattice also removing itinerant 
charge carriers, and (iii) application of both Zn substitution and excess oxygen results in a com- 
pensation effect as Zn substitution partially removes the Fe vacancies. Both electrical conduction 
and magnetism is determined by the density and hopping amplitude of the itinerant charge carriers 
on the B sublattice, providing electrical conduction and ferromagnetic double exchange between 
the mixed- valent Fe^/Fe^ - ions on the B sublattice. A decrease (increase) of charge carrier den- 
sity results in a weakening (strengthening) of double exchange and thereby a decrease (increase) of 
conductivity and the saturation magnetization. This scenario is confirmed by the observation that 
the saturation magnetization scales with the longitudinal conductivity. The combination of tailored 
Zn a: Fe3_ 2 ,04 films with semiconductor materials such as ZnO in multi-functional heterostructures 
seems to be particularly appealing. 

PACS numbers: 75.70.-i 81.15.Fg, 85.75.-d 75.50.Dd 



1. INTRODUCTION 

Spintronic materials and devices are in the focus of 
current research activities 1,2 . Regarding functional spin- 
tronic materials, transition metal oxides are promising 
candidates, since they show a rich variety of electrical and 
magnetic properties specifically interesting for applica- 
tions in spintronics. Evidently, useful spintronic devices 
such as magnetoresistivc elements based on the tunnel- 
ing magnetoresistance 3 or spin transistors 4 require fer- 
romagnetic materials with high Curie temperature Tc 
well above room temperature and large spin polariza- 
tion P of the charge carriers at the Fermi level E-p. 
Moreover, advanced applications often require materi- 
als with magnetic properties which can be deliberately 
tuned by external control parameters such as an applied 
electric field 5 , elastic stress 6 ' 7 , or light 8 . In this respect, 
Fc304 (magnetite) is a highly promising candidate. First, 
it has both a high Curie temperature of Tq ~ 860 K 
and according to band structure calculations is expected 
a half-metal 9 corroborated by spin-resolved photoelec- 
tron spectroscopy 10,11 . Second, recent experiments indi- 
cated that both the electronic and magnetic properties 
of FC3O4 thin films can be nicely tailored in solid solu- 



tion systems of Fc 3 - x M x 4 (M = Zn, Mn) 12 ' 13 . Third, 
it has been demonstrated that the magnetic anisotropy 
of Fe304 can be tuned by clastic stress imposed by a 
piezoelectric actuator 14 . Taken all together, due to their 
interesting, versatile, and tunable properties magnetite 
thin films and heterostructures became the focus of re- 
cent research activities. 

Before discussing the various possibilities for tuning 
and tailoring the electronic and magnetic properties of 
magnetite we briefly review its basic structural and mag- 
netic properties. Magnetite is known to have an inverse 
spinel structure as shown in Fig. 1(a). The A sites (8 
per unit cell), which are surrounded by oxygen tetrahc- 
dra, are occupied by trivalent Fe 3 ^ ions (3d 5 , S = 5/2), 
whereas on the octahcdrally coordinated B sites (16 per 
unit cell) there is an alternating arrangement of Fe 2 ^ 
(3d 6 ,S = 2) and Fe| + (3d 5 ,S = 5/2) ions. There- 
fore, the sum formula of magnetite can be expressed as 
[Fe 3+ ] J 4[Fe 3+ Fe 2+ ]_B04. The density of itinerant charge 
carriers is determined by the density of the tig spin-down 
electron on the B site, i.e. by the density of Fe 2 ^. The 
magnetic exchange in magnetite is governed by a combi- 
nation of antiferromagnetic superexchange (SE) and fer- 
romagnetic double exchange (DE) interactions. There 
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are three antiferromagnetic SE interactions Jaa (A-O- 
A), J BB (B-O-B), and J AB (A-O-B) between the Fe 3+ 
ions on the A and B sites mediated by the oxygen (O) 
ions. In addition, there is a ferromagnetic DE interac- 
tion mediated by the itinerant spin-down ti g electrons 
hopping between the mixed-valent Fe ions on the B sites 
(cf. Fig. 1). Owing to Hund's rule coupling, the spins 
of these itinerant electrons are antifcrromagnctically cou- 
pled to the localized spins formed by the 3d spin-up elec- 
trons. As pointed out by Neel 15 , in the simplest model 
fcrrimagnctism in Fe 3 04 with high Tq is obtained with- 
out any DE interaction for large Jab 3> Jbb, Jaa forcing 
an antiparallel alignment of the moments on the A and B 
sites (see Fig. 1(b)). Since the antiparallel Fe 3 A + and Fe^ + 
moments compensate each other, a saturation magneti- 
zation of 4^iB/formula unit (f.u.) is expected from the 
remaining Fc 2 B + (S = 2) moments. Later on, the sim- 
ple Neel model has been extended by Yafet and Kittel 16 . 
They proposed a more elaborate model in which the B 
sublattice is subdivided into two Fe^ + and Fe^ + sublat- 
tices. It was shown that on weakening Jab and strength- 
ening Jbb , the B site magnetic moments arc no longer 
rigidly parallel to the A site moments. The stronger B- 
O-B SE interaction results in spin canting expressed by a 
finite Yafct-Kittcl angle and thus a reduction of the sat- 
uration magnetization. More recent models show that 
a detailed modelling of the magnetic properties of mag- 
netite is only possible by taking into account the ferro- 
magnetic DE interaction of the B sublattice competing 
with the antiferromagnetic SE 17,18,19 . 

There are several possibilities to tailor the magnetic 
properties of magnetite. First, Fe A + on the A site can 
be replaced by an isovalent nonmagnetic ion. In this way 
antifcrromagnctically coupled moments on the A sites are 
removed without affecting the magnetic exchange on the 
B sublattices. Hence, an increase of the saturation mag- 
netization is expected at low substitution level. How- 
ever, on increasing the A site substitution the dilution 
of the A site moments reduces Jab, whereas Jbb stays 
about constant. Then the moments on the A and B 
site do no longer need to be rigidly antiparallel, result- 
ing in a finite Yafet-Kittel angle (see Fig. 1(c) and (d)). 
This angle increases with increasing substitution, result- 
ing in a reduction of the saturation magnetization. Due 
to the isovalent substitution the charge carrier density 
should stay unaffected. However, the hopping between 
the mixed-valent (Fe 2 B + /Fe 3 ^) ions on the B sublattice is 
reduced due to the spin canting, resulting in a lower elec- 
trical conductivity and reduced DE interaction. To our 
knowledge a nonmagnetic isovalent A site substitution in 
Fe3C>4 has not been reported so far. Most likely this can 
be attributed to the fact that suitable trivalent ions usu- 
ally prefer the octahedral coordination on the B site and 
thus cannot be substituted solely on the tetrahedral A 
site. Second, as done in our work, Fe 3 A + on the A site can 
be replaced by a nonmagnetic divalent ion such as Zn 2+ 
(3d 10 , 5 = 0) 20 . It is known that Zn 2+ preferably occu- 
pies the tetrahedrally coordinated A site in the inverse 
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FIG. 1: (color online) (a) Inverse spinel structure of mag- 
netite with Fe^ + on the tetrahedrally coordinated A site and 
mixed-valent Fe|+ (3d 6 , S=2)/Fe|+ (3d 5 ,S = 5/2) ions on 
the octahedrally coordinated B site. Zn substitution replaces 
Fe A + on the A site by nonmagnetic Zn 2+ (3d 10 , 5 = 0) ions. 
Also shown is the occupation of the Fe^ + and Fe^ + 3d states 
with t2 g and e 9 symmetry separated by the crystal field split- 
ting due to the octahedral oxygen surrounding. The t2 g spin- 
down electron can hop between Fe^ + and Fe^ + ions with its 
spin coupled anti-parallel to the local moment formed by the 
spin-up electrons owing to Hund's rule coupling, (b) - (e) 
shows the A and B sublattice configuration for Zn a; Fe3- 2 ,04 
for (b) x = 0, (c) x = 1/8, (d) x = 0.5, and (e) x = 1. 



spinel structure 21 . Again some of the antiferromagncti- 
cally coupled moments on the A sites are removed, lead- 
ing to an initial increase of the saturation magnetization 
at low Zn substitution 12,22 . However, in the same way as 
discussed before, a finite Yafet-Kittel angle and, in turn, 
a reduction of the saturation magnetization is expected 
going to larger substitution levels 22 . Moreover, substi- 
tution of Fe 3 A + by Zn 2+ on the A site also reduces the 
amount of Fe 2 B + on the B site due to charge neutrality. 
That is, the amount of itinerant charge carriers mediat- 
ing the DE on the B sublattice is reduced as shown by 
photoemission spectroscopy 13 . Thus, with increasing Zn 
substitution the electrical conductivity is reduced both 
by a reduction of the density of itinerant charge carriers 
and their hopping amplitude due to spin canting. Third, 
Fe vacancies can be introduced by preparing magnetite 
samples in excess oxygen 20,23 ' 24 ' 25,26 . Since Fe vacan- 
cies are formally equivalent to the presence of excess of 
2 ~ ions, charge neutrality again requires a reduced (in- 
creased) amount of Fe 2 ^ (Fe 3 ^) ions on the B site. This 
strengthens antiferromagnetic SE and weakens ferromag- 
netic DE on the B sublattice. In turn, this results in a 
finite Yafct-Kittcl angle and a reduced saturation magne- 
tization. Furthermore, a reduced electrical conductivity 
is expected both by a reduced carrier density and hop- 
ping amplitude. 

Our brief discussion shows that Fe based magnetic 
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FIG. 2: (color online) X-ray diffraction diagrams (lj-29 scans) of the out-of-plane reflections from Zn a; Fe3_ a; 04 thin films 
deposited in different growth atmospheres of Ar/02 (99:1) (left panels) and pure Ar (right panels). With increasing x the 
Zn x Fe3_ a: 04 (004) reflection shifts across the (002) reflection of the MgO substrates marked by the vertical line. The insets 
show the corresponding rocking curves of the Zn^Fea-^C^ (004) reflection, except for (f) where the Zn^Fea-^CU (008) reflection 
is displayed. 



oxides with spinel structure are interesting materials. 
They do not only have high Curie temperature and 
high spin polarization but also offer various opportu- 
nities to tailor their electrical (charge carrier density, 
electrical conductivity) and magnetic properties (satu- 
ration magnetization, Curie temperature). Such tun- 
able ferromagnetic materials are desired for spintronic 
devices operating at room temperature. For example, 
tuning the electrical conductivity of ferromagnetic ma- 
terials with high spin polarization is promising for spin 
injection devices as one can reduce and deliberately 
choose the conductivity mismatch between semiconduc- 
tors and ferromagnetic spin injectors. To this end, re- 
cently the growth of epitaxial FcaCj/ZnO heterostruc- 



tures has been demonstrated 27 . Here, we report on the 
growth as well as the structural, magnetic and magne- 
totransport properties of epitaxial thin film samples of 
Zn :r Fe 3 _ ;l .04 with x = 0,0.1,0.33,0.5,0.9 deposited in 
different oxygen partial pressure. We show that their 
saturation magnetization and electrical conductivity can 
be tuned over a wide range both by Zn substitution and 
varying oxygen partial pressure during growth. Further- 
more, our systematic study provides clear evidence that 
their electrical conductivity is closely related to the over- 
all magnetization. We also carefully consider the pres- 
ence of anti-phase boundaries (APBs) in the thin film 
samples which could affect both their electrical conduc- 
tivity and their saturation magnetization 28 ' 29 . However, 
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we have clear evidence that APBs do not play a dom- 
inant role in this study. In particular, they cannot ex- 
plain the change of the magnetic and magnetotransport 
properties on varying the Zn content and the deposition 
atmosphere. In contrast, these properties can be consis- 
tently explained solely by disorder and spin canting on 
the B sublattice, tending to localize the itinerant charge 
carriers on the B sublattice. 



2. THIN FILM GROWTH 

Epitaxial thin films of Zn a ,Fe3_ a; 04 with x < 0.9 and 
thicknesses between 40 nm and 60 nm were deposited 
from stoichiometric targets on single crystalline, (001) 
oriented MgO substrates by laser molecular beam epi- 
taxy (laser-MBE) 30 . The exact film thickness was deter- 
mined by x-ray reflectometry. Both bulk Fe304 having an 
inverse spinel structure (Fd3m) and MgO {Fmim) are 
cubic with lattice constants of aFe 3 o 4 

= 8.394 A 31 and 

fl MgO = 4.212 A, respectively, resulting in a small lattice 
mismatch of (aFe 3 o 4 — 2aM g o)/2aMgO = —0.4%. The en- 
ergy density of the KrF excimer laser (A = 248 nm) at the 
target was 3.1J/cm 2 and the laser repetition rate 2 Hz. 
At the same total pressure of 3.7 x 10~ 3 mbar, two sets 
of samples were deposited in two different growth atmo- 
spheres. In pure Ar, thin films with x = 0, 0.1, 0.5, and 
0.9 were grown at a substrate temperature of 320°C. In 
an Ar/02 (99:1) mixture, thin films with x = 0.5 were de- 
posited at 320°C and with x = 0, 0.33, and 0.9 at 400°C. 
The growth process was monitored in-situ by reflection 
high energy electron diffraction (RHEED) 32 . We ob- 
served four RHEED intensity oscillations per unit cell, in- 
dicating a block-by-block growth mode with four charge 
neutral blocks. This is already known for Fe3 4 33,34 and 
also holds when substituting Zn up to Zn .gFe 2 .iO4. A 
more detailed description of the thin film deposition for 
Fe304 on MgO is given elsewhere 33 ' 35 . 



3. STRUCTURAL PROPERTIES 

The structural properties of the Zti\ x ¥c^^ x Oi films 
were analyzed by high resolution x-ray diffractometry us- 
ing a Bruker AXS D8 Discover four-circle diffractometer. 
Our detailed analysis reveals a very high epitaxial qual- 
ity of all thin film samples. In [00^] direction, ui-29 scans 
display no impurity phases - in particular from other 
iron or zinc oxides. We note that we may not be able 
to distinguish between Fe304 and 7-Fe203 (maghemite) 
as they both share the same inverse spinel lattice 36 . It 
is known that excess oxygen results in the formation 
of Fe vacancies and in the extreme case to formation 
of Fe 3 04/7-Fe 2 03 solid solutions 24 ' 25 . The Fe 3 4 /7- 
Fe203 thermodynamic equilibrium line was determined 
to log( Po2 /Po) = - 24 'j 4K with p Q = 9.1 x 10 13 bar 24 . 
Unfortunately, both calculations and experimental stud- 
ies of phase and point defect equilibria always apply to 
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FIG. 3: (color online) (a) Change of the c axis lattice pa- 
rameter of Zn a; Fe3_ a ;04 films grown on MgO(OOl) in pure Ar 
atmosphere (squares) and an Ar/02 mixture (circles) with 
Zn content x. The bulk c axis lattice constant of Fe304 is 
marked by the dashed horizontal line, (b) Reciprocal space 
map around the (204) reflection of the MgO(OOl) substrate. 
The (408) reflection from the Zno.1Fe2.9O4 film appears at 
the same in-plane scattering vector marked by the dashed 
vertical line. 



thermal equilibrium situations which certainly are not 
appropriate for the laser-MBE growth process occuring 
far from equilibrium. Nevertheless, the high energy of the 
particles in the laser plume can be associated to an effec- 
tive growth temperature much larger than the substrate 
temperature favoring the stability of Fe3 04 24,25 . There- 
fore, for the deposition temperature and oxygen partial 
pressure used in our experiments we are sufficiently far 
away from the Fe304/7-Fc203 equilibrium line 37 - 38 . 

As shown by Fig. 2, the (004) reflection from 
Zn 2; Fe3_ 2 ;04 can be clearly observed together with satel- 
lites due to Laue oscillations. These oscillations demon- 
strate that the thin films are coherently strained and 
have a small surface roughness. The rocking curves of 
the (004) or (008) reflections (insets of Fig. 2) show a 
full width at half maximum (FWHM) of Auj < 0.05°. 
This demonstrates the excellent structural quality of 
the Znj;Fc3_ x 04 thin film samples with very low mo- 
saic spread, comparable to epitaxial FC3O4 films 33,35 . In 
the ui-28 scan, the (004) reflection from Fe304 is posi- 
tioned at a slightly larger angle than the (002) reflec- 
tion from the MgO substrate (see Fig. 2(a,b)), indicating 
CFe 3 4 < 2cM g o- With increasing Zn concentration, the 
(004) reflection moves to lower angles and crosses the 
MgO(002) reflection at x = 0.5 (see Fig. 2(e,f)). At 
x = 0.9, the position of the (004) reflection from the 
Zn a; Fe3_ 2 ;04 film is shifted to an angle below the (002) 
reflection from MgO (see Fig. 2(g,h)). This corresponds 
to CZno.9F02.1O4 > 2cM g o- 

The absolute value and variation of the c axis lattice 
parameter of the Zn a; Fe3_ a; 04 films with Zn content x are 
related to two effects, namely epitaxial coherency strain 
and the larger ionic radius of tetrahcdrally coordinated 
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Zn 2+ of 0.6 A compared to the radius of only 0.49 A of 
Fe 3+ on the tetrahedrally coordinated A site 12 . First, for 
x = the lattice mismatch of —0.4% between Fe3C>4 film 
and MgO substrate leads to a tensile in-plane epitaxial 
coherency strain and, in turn, to a slight reduction of 
the c axis lattice parameter of the Fe3C>4 film below the 
bulk value (cf. Fig. 3(a)). The fully coherent growth of 
the Zn x Fe3- x Oi thin films is demonstrated by recipro- 
cal space maps around the (204) reflection of the MgO 
substrate (see Fig. 3(b)). The (408) reflection of the film 
is located exactly at the same in-plane scattering vector 
qn as the (204) reflection of the substrate. This clearly 
demonstrates that the in-plane lattice constant of the film 
perfectly matches twice the one of the MgO substrate. 
Second, for x > the unit cell volume of Zn r Fe3_ ;E 04 
is expected to increase about linearly with increasing x 
due to the larger ionic radius of Zn 2+ . However, since the 
Zn K Fe3_2;04 films grow coherently on the MgO(OOl) sub- 
strate without any relaxation up to the maximal thick- 
ness of 60 nm studied in our experiments, the in-plane 
lattice constant stays unchanged. Only the out-of-planc 
lattice constant is found to vary with increasing x. De- 
pending on whether the in-plane strain is tensile or com- 
pressive, a reduction or expansion of the unit cell in [00£] 
direction, respectively, is found (Fig. 3(a)), resulting in 
a tetragonal distortion of the cubic lattice. The c axis 
parameter increases from a value below {x = 0, tensile 
strain) to a value above the bulk value {x > 0.3, com- 
pressive strain). The expansion of the out-of-plane lattice 
parameter with increasing x has been reported also for 
(lll)-oriented films grown on Al 2 O3(0001) substrates 12 . 
However, we note that those films are relaxed due to 
the very large lattice mismatch of 8% between film and 
substrate. Therefore, both results cannot be compared 
directly. 

An interesting result shown in Fig. 3 is the fact that 
the c axis lattice parameters of the Zn a: Fe3_ 2 ;04 films 
grown in pure Ar atmosphere are larger than those of 
the films grown in an Ar/02 mixture and also show a 
slightly weaker increase with Zn content in the range 
< x < 0.5. This most likely is caused by the formation 
of Fe vacancies when growing the films in finite oxygen 
partial pressure 20 ' 23,24 ' 25,26 . The presence of Fe vacancies 
can be viewed as an internal negative pressure effect lead- 
ing to a reduced cell volume. The weaker increase of the 
c axis parameter for the films grown in Ar/02 suggests 
a partial compensation of the Fe vacancies by Zn sub- 
stitution. This is in agreement with the transport data 
discussed below. We note, however, that this result has 
to be considered with some care. The problem is that 
the close vicinity of film and substrate reflections in the 
lj-29 scans does not allow an unambiguous derivation of 
the film c axis parameter with small error bars. 

In summary, the structural analysis demonstrates that 
our Zn a; Fe3„ a ;O4(001) thin films are coherently strained 
and show a very small mosaic spread. Zn substitution 
results in an increase of the unit cell volume, causing the 
in-plane epitaxial strain to change from tensile to com- 



pressive on increasing the Zn content x. Growth in excess 
oxygen results in a reduction of the unit cell volume most 
likely due to the formation of Fe vacancies. The latter 
most likely are partially compensated by additional Zn 
substitution. 



4. MAGNETOTRANSPORT PROPERTIES 

For magnetite the electrical conduction above the Ver- 
wey transition is believed to be determined by the hop- 
ping of the spin-down electrons between the mixed- valent 
Feg and Fc 3 + ions on the B sublattice. From the hop- 
ping amplitude t ~ 0.1 eV an electron conduction band- 
width D ~ 1 eV is expected. However, experimental data 
cannot be explained within a simple band theory. They 
can be better ascribed to some thermally activated mo- 
tion of charge carriers. The reason is that the bare hop- 
ping amplitude and bandwidth is strongly reduced due 
to the translation of a polaronic lattice deformation as- 
sociated with electron motion 39,40 . If the reduced band- 
width is smaller than the energy of polarization phonons 
(typically 70meV in magnetite), small polaron hopping 
is expected 41 . We note that the hopping takes place be- 
tween Fe 2 ^ and Fe 3 ^ ions on the B sublattice. Below 
the Curie temperature, the spins of these ions arc aligned 
parallel in the ideal case. However, as discussed above in 
reality there is some spin canting on the B sublattice. 
This has to be taken into account by an appropriate spin 
correlation factor 42 leading to a spin dependent part of 
the activation energy. Since the spin alignment is im- 
proved by applying a magnetic field, this spin-dependent 
part of the activation energy can be reduced by the exter- 
nal field resulting in a large negative magnetoresistance 
as shown below. 

For magnetotransport measurements, the films were 
patterned into typically 45 iim wide and 350 /im long Hall 
bars directed in the (010) direction by photolithography 
and Ar ion beam milling 43 . The longitudinal resistivity 
p xx has been measured as a function of temperature T 
and applied magnetic field H using a standard four-probe 
technique. The films with x = 0.9 were found to be 
insulating and are not further discussed in the following. 
We note that measurements of the anomalous and the 
ordinary Hall effect have been reported previously 35,44 
and arc not discussed here. 



4.1. Temperature dependence of resistivity 

The p xx (T) curves of representative samples with x < 
0.5 are shown in Fig. 4. Despite a similar shape of all 
Pxx(T) curves, it is evident that the absolute magni- 
tude of p xx sensitively depends on the Zn substitution 
x and the growth atmosphere. For both samples sets 
grown in pure Ar atmosphere and an Ar/02 (99:1) mix- 
ture, p xx (T) increases by more than two orders of magni- 
tude on decreasing T from 375 K down to 90 K. Plotting 
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\og{p X x) on a reciprocal T scale as shown in Fig. 4(b) 
results in about linear curves over a wide temperature 
range except for the x = and x — 0.1 samples grown 
in pure Ar atmosphere. The kink in p xx (T) of the x = 
sample results from the Verwey transition discussed be- 
low. Fig. 4(b) suggests that the longitudinal resistivity of 
the samples with x > 0.1 follows a simple activated be- 
havior Pxx(T) oc exp(£? a /£;BT) with the Boltzmann con- 
stant k-Q = 1.38 x 10 -23 J/K and activation energies E a 
ranging between 61 and 84meV. However, an unambigu- 
ous determination of the detailed transport process is dif- 
ficult. As indicated by Fig. 4(c), the p X x{T) data also can 
be equally well fitted by small polaron hopping 41,45,46 , 
Pxx{T) oc T exp{E vo \/kBT), with the potential barrier for 
polaron hopping, E po \, ranging between 64 and 95meV. 
Actually, small polaron hopping yields a slightly bet- 
ter fit for the samples grown in pure Ar atmosphere, 
whereas the simple activated transport better fits the 
transport data of the samples grown in an Ar/02 mix- 
ture. The derived activation energies are typical for zinc 
ferrites 20,37,46,47,48 . We also note that a similar transport 
behavior with similar activation energies is reported for 
the mixed-valent manganitcs, where charge transport is 
determined by the hopping between mixed-valent Mn 3+ 
and Mn 4 + i on s 49,50,51,52 . 

We next discuss the observability of the Verwey 53 tran- 
sition in the measured p xx {T) curves. In bulk mate- 
rial, this metal-insulator transition, associated with a 
cubic to monoclinic structural phase transition, occurs 
at Ty = 123 K and in the classic picture is argued 
to arise from a charge ordering of the Fe 2 ^ and Fe 3 ^ 
on the B sublattice in a process reminiscent to Wigner 
crystallization 54 . Whereas at T > Ty the extra electrons 
from Fe 2 ^ ions can hop to neighboring Fe 3 ^ ions on the 
corner sharing tctrahedral network on the B sublattice, 
giving rise to electrical conduction, this process is frozen 
in at T < Ty due to Coulomb repulsion 55 . However, 
more recently it became evident that also elastic and or- 
bital interactions play a significant role. This leads to a 
renewed controversial discussion of the nature of the Ver- 
wey transition 18,56,57,58,59 . It is evident from Fig. 4 that 
a pronounced change in the p xx (T) curves due to the Ver- 
wey transition is observed only for the Fe3C>4 film (x = 0) 
grown in pure Ar atmosphere. For this film there is no 
jump but a clear kink in the p X x(T) curve at a slightly 
reduced temperature of about 115 K. For all other sam- 
ples no clear signature of the Verwey transition is ob- 
servable. This can be understood in a straightforward 
way. As discussed above, both substitution of Zn and 
introducing Fe vacancies by growing the films in finite 
oxygen partial pressure results in deviations from equal 
numbers of Fe 2 ^ and Fe 3 ^ on the B sublattice. This 
can be viewed as disorder on the B sublattice, tending to 
smear out the Verwey transition. The fact that the ob- 
servability of the Verwey transition seems to sensitively 
depend on deviations from perfect stoichiomctry already 
has been reported in literature 20,23,25,26,28 and is further 
supported by our results. The other way round, we can 
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FIG. 4: (color online) (a) Longitudinal resistivity p xx as 
a function of temperature T for epitaxial Zn I Fe3_ a: 04 thin 
films, grown in pure Ar atmosphere (solid lines) or an Ar/02 
(99:1) mixture (dashed lines) on MgO(OOl) substrates. In 
(b) and (c) p xx and p X x/T of the same samples are plot- 
ted versus a reciprocal temperature scale. Fitting the linear 
parts of the curves in (b) and (c) by p xx oc exp(E a /kBT) and 
p xx /T oc exp(_Z5p i/fcsr), respectively, gives the listed activa- 
tion energies E a and E po \. 



argue that the absence of any clear signature of the Ver- 
wey transition for the Fe3C>4 films substituted by Zn or 
grown in finite oxygen partial pressure provides evidence 
for the presence of an unequal number of Fe 2 ^ 
ions on the B sublattice. 



2+ and Fe|+ 



We next discuss the p xx values of the Zn a; Fe3_ 2 ,04 
(0 < x < 0.5) films grown in pure Ar atmosphere. The 
strong increase of resistivity with decreasing T can be un- 
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derstood in terms of a thermally activated hopping type 
transport mechanism of the itinerant electrons on the B 
sublattice. For x = 0, the room temperature p xx value 
(7.6 x 10~ 5 57m) corresponds well with literature data 
for Fe3C>4 thin films 35 or single crystals 60 . In this con- 
text, we note that the presence of anti-phase boundaries 
(APBs) in Fe3C>4 thin film samples may have a signifi- 
cant influence on the measured longitudinal resistivity 61 . 
However, this is not the case for our samples. First, our 
Fe304 thin film grown in Ar atmosphere has a room- 
temperature resistivity value close to that of Fe3C>4 single 
crystals (6.2 x 10~ 5 flm) considered as APB-free 60 . Sec- 
ond, comparing the resistivity value to those of Fe3C>4 
thin films with different APB density, we estimate that 
the volume fraction of the APB phase in our thin film 
is less than 20% 61 . Third, it is known that the forma- 
tion and the density of APBs sensitively depends on the 
used substrates 62 or the film thickness 61 . However, in 
our study always the same substrate and film thickness 
is used. Since also Zn substitution is not expected to 
severely affect the APB density, we can safely assume 
that all samples grown in the same atmosphere have a 
similar APB density. 

With increasing x, as discussed in the introduction, 
Zn 2+ preferably occupies the tetrahedrally coordinated 
A site in the inverse spinel structure 21 . Hence, the sub- 
stitution of Fe A + by Zn 2+ reduces the amount of Fe 2 ^ 
as 

[Fe 3 +UFe 3 +Fc 2 +] B 2 -^ 

[Zn^+FeflJ^ [Fe^Fc^^O 2 -. (1) 

That is, the density of itinerant electrons on the B sublat- 
tice is reduced, resulting in an increase of the resistivity 
with increasing x. This is in perfect agreement with our 
observation and literature data 12,13 . We note however, 
that the increase of resistivity is not only caused by a 
reduction of the carrier density but also by several other 
effects. First, deviations from the 1:1 balance of Fe 2 ^ 
and Fe 3 ^ on the B sublattice causes disorder, tending 
to localize the itinerant charge carriers. Second, Zn sub- 
stitution dilutes the A sublattice, thereby weakening the 
antifcrromagnctic exchange Jab between the A and B 
sublattices. As discussed above, this results in spin cant- 
ing on the B sublattice, reducing the hopping amplitude. 
This amplitude is maximum for a parallel alignment of 
the local moments due to strong antiparallel Hund's rule 
coupling of the spins of the itinerant electrons. Both 
effects arc expected to result in an increase of the acti- 
vation energy for the hopping transport with increasing 
x in agreement with our observations. 

The p xx (T) curves of the Zn 2 ,Fe 3 _ 2; 04 (0 < x < 0.5) 
films grown in an Ar/02 mixture show a similar overall 
temperature dependence, however, with resistivity val- 
ues that are much higher than those measured for the 
films grown in pure Ar atmosphere. Furthermore, the 
resistivity values decrease with increasing x in contrast 
to what is observed for the films grown in Ar. We first 



discuss the origin of the much higher resistivity value 
( Pxx = 7 x 1(T 3 Om at 300 K) of the Fe 3 4 (x = 0) film 
grown in finite oxygen partial pressure compared to that 
(p xx = 7.6 x 10~ 5 ftm at 300 K) of the film grown in 
pure Ar. Since growth of magnetite in excess oxygen is 
known to create Fe vacancies 

X 023,24,25,26 j equ i va l cnt to 

two (three) missing electrons per Fe 2+ (Fe 3+ ) vacancy, 
the requirement of charge neutrality shifts the 1:1 bal- 
ance between Fc 2 ^ and Fc 3 ^ towards Fe 3 ^ as 
[Fe 3 +UFe 3+ Fc?+] B 2 -^ 

[Fe? + ] A [Fe 3 + 2 ,Fe 2 l 35 X°] B 2 -. (2) 

Evidently, Fe vacancies in the same way as Zn substi- 
tution result in a reduction of Fc 2 ^. Hence, the den- 
sity of itinerant electrons on the B sublattice decreases 
and, in turn, the resistivity increases with increasing 6. 
This is exactly what is observed in our experiments and 
reported in literature 20 ' 63 . Actually, for Zn a; Fe3_ 2 ,04 
(0 < x < 0.3) single crystals Wang et al. report a compo- 
sitional correspondence x ^ 35 regarding the measured 
resistivity data 20 . Again, with the same arguments given 
above for the films grown in pure Ar, a reliable evaluation 
of the various mechanisms (localization, spin canting) re- 
sponsible for the increase in resistivity is not possible. We 
also note that according to eq.(2) growth in high oxygen 
partial pressure (not the case in our experiments) may 
lead to the formation of the cubic 7-Fe2C>3 (maghemite) 
phase in the extreme case of 5 = 1/3 64 . 

Zn substitution in our magnetite films grown in an 
Ar/02 mixture results in a decrease of resistivity in con- 
trast to what is observed for the films grown in pure Ar 
atmosphere. On first sight, this is astonishing and seems 
to be in conflict with the above discussion, since Zn 2+ 
again preferably occupies the tetrahedrally coordinated 
A lattice sites 21 , substituting Fe 3+ and thereby reducing 
carrier density. That is, one would expect an increase of 
resistivity with increasing x. However, since Zn substi- 
tution also results in an increase of the unit cell volume 
due to the larger ionic radius of Zn 2+ compared to Fe 3+ , 
it is expected that Zn substitution removes part of the 
Fe vacancies. This is intuitive since now Fe can be more 
easily incorporated into the expanded lattice. Then, it is 
expected that Zn 2+ does not only remove an electron by 
substituting Fe 3 A + on an A site, but also adds two/three 
electrons by removing a vacancy on a Fe^/Fe 3 ^ site. 
That is, by Zn substitution also X° is replaced by Fe 2+ or 
Fe 3+ . Evidently, in total this leads to an effective increase 
of the carrier density with increasing x in agreement with 
our experimental data in the range up to x = 0.5. The 
partial removal of the Fe vacancies by Zn substitution is 
in agreement with the x-ray data of Fig. 3(a). 

In summary, the p xx (T) curves of the Zn a; Fe3_ :r 04 
films grown under different oxygen partial pressure show 
a similar shape, originating from a hopping type trans- 
port mechanism, but strongly differing absolute resis- 
tivity values. These differences can be consistently ex- 
plained by the change of the itinerant charge carrier den- 
sity on the B sublattice and their hopping amplitude by 



8 



either Zn substitution or the generation of Fe vacancies 
due to growth at finite oxygen partial pressure. How- 
ever, applying both Zn substitution and excess oxygen 
results in an increase of the carrier density as Zn substi- 
tution is removing part of the Fe vacancies generated by 
excess oxygen. It would be interesting to directly check 
the change of the carrier density by measurements of the 
Hall effect. However, we cannot unambiguously separate 
the small ordinary from the large anomalous Hall contri- 
bution and any attempt to do so would result in large 
errors of the derived carrier density 44 . 

We note that the observed dependence of the resistiv- 
ities on the Zn substitution levels x cannot be simply 
explained by the assumption of different APB densities 
in the different thin film samples. It is unlikely that with 
increasing x the APB density increases for samples of 
the same thickness when grown in Ar atmosphere while 
it decreases for those prepared in the Ar/02 mixture. 
Moreover, the samples show a universal scaling relation 
of the anomalous Hall conductivity with the longitudinal 
conductivity indicating a negligibly small APB resistivity 
as already pointed out earlier 44 . 
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4.2. Magnetoresistance 



In this subsection we address the magnetoresistance 
(MR) of the Zn^Fca.^CU films. The MR effect, MR la . = 
[Pxx(H) - Pxx(0)]/Pxx(0) is shown in Fig. 5(a) at several 
temperatures for the x = 0.1 samples grown in pure Ar. 
As shown in Fig. 5(b), very similar MR XX (H) curves with 
similar absolute values of MR II; at room temperature are 
found for the other samples, although the absolute val- 
ues of their resistivities vary by more than two orders of 
magnitude (cf. Fig. 4). This interesting observation can 
be consistently explained in the framework of thermally 
activated hopping of the itinerant electrons on the B sub- 
lattice. Due to the strong on-site Hund's rule coupling, 
the spins of the itinerant spin-down electrons are coupled 
antiparallel to the localized spins formed by the spin-up 
electrons. Therefore, the activation energy for the hop- 
ping process is given by the sum of a spin independent 
energy (E a or E po i) and a spin dependent contribution 
E s . The latter depends on the spin canting of the local 
moments on the B sublattice and disappears for a per- 
fect parallel alignment. Since in the presence of a finite 
canting of the local moments an applied magnetic field 
tends to improve the alignment, a magnetic field depen- 
dent total activation energy E(H) = E a , po i+E s (H) is ob- 
tained. Then, no matter whether the transport is by ther- 
mally activated hopping, p xx = pQ^ a exp[E(H)/kBT] 1 or 
small polaron hopping, p xx /T = p 0tPO icxp[E(H)/k B T] 7 



FIG. 5: (color online) (a) MR„ plotted versus the magnetic 
field H applied perpendicular to the film plane at different 
temperatures for Zno.1Fe2.9O4 thin films grown in pure Ar 
atmosphere, (b) MRn (H) curves measured at 300 K for dif- 
ferent Zn^Fes-zCU grown in pure Ar atmosphere (left) or an 
Ar/02 mixture (right). The dash-dotted gray lines are fits to 
the data according to eq.(6). 



the MR effect is obtained as 



MR XI = 




with AE S (H) = E S (Q) - E S (H). Evidently, with the as- 
sumptions made, MR XX (H) is independent of the magni- 
tude of the prefactors (/?o, a or po.poi) and the activation 
energies {E a or -E po i), which may strongly vary for sam- 
ples with different x and grown in different atmosphere. 
Moreover, similar MR XX (H) curves with comparable ab- 
solute values are obtained, if the magnetic field induced 
change AE S (H) of the spin-dependent part of the acti- 
vation energy is about the same for all samples. 

Since the transport in Zn :E Fe3_ a ;04 by hopping of the 
t2 g spin-down electrons between mixed- valent Fe^ + /Fe^ + 
ions is equivalent to the hopping of the e g electrons be- 
tween mixed- valent Mn 3+ /Mn 4+ ions in the doped man- 
ganites, we can adopt models developed for the MR effect 
in the perovskite manganitcs 65,66 ' 67 . According to these 
models the magnetic field dependent change of the acti- 
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vation energy for hopping between two lattice sites i and 
j can be expressed as 

AE S (H) = a ((mi • m.j)H - (mi • mj}jj=o) , (4) 

where m, and nij are the local moments on the lattice 
sites i and j, (. . .)h=o and (■ ■ -)h denote the sample av- 
erage at zero and finite magnetic fields, respectively, and 
a is a constant that may depend on temperature. The 
local moments m^ and m.j can be associated either with 
the localized spins of the individual Fe ions or the mo- 
ments of small spin clusters. Note that these spin clusters 
also can be considered as spin polarons, that is, itinerant 
electrons dressed by a small cloud of parallel aligned lo- 
calized spins. Evidently, AE S = is independent of H, 
if the local moments are already perfectly aligned in zero 
magnetic field (ideal fcrromagnct). In this case, no MR 
effect is expected. However, if there is a finite canting 
of the neighboring moments, applying a magnetic field 
results in a reduction of this spin canting and thereby 
in AE S > 0. According to cq.(3) this results in a finite 
negative MR effect with a magnitude determined by the 
value of AEs/ksT. The magnetic field dependence of 
AE S can be estimated by keeping in mind that in the 
ferromagnetic state both the molecular field and the ap- 
plied magnetic field support a parallel alignment of the 
local moments. Since the molecular field is usually much 
larger than the applied field, the deviations of the di- 
rections of the local moments from the average direction 
are small. In this case AE S (H) is found to follow the 
Brillouin function 03 giving 65 ' 66 ' 67 

AE S (H) = /J*(£g2?) . (5) 

with the vacuum permeability /j.q = 4ir x 10~ 7 Vs/Am. 
Here, // e ff is the average value of the local moments and 
(3 a constant, which may show a weak temperature de- 
pendence. Assuming AE S <C fc^T, what certainly is a 
reasonable assumption at room temperature, we can ap- 
proximate eq.(3) by 



MR, 




This shows that the magnetic field dependence of the 
MR effect measured at constant temperature should fol- 
low the Brillouin function. As shown in Fig. 5(b), this 
is indeed the case. The data of all samples can be well 
fitted by eq.(6). Fitting the data gives the n e s values 
of the local moments. They range between 70 to 120 /ib 
and 60 to 70 /is, corresponding to about 18 to 30 and 15 
to 18 Fe 2+ ions, for the samples grown in pure Ar and 
an Ar/C>2 mixture, respectively. Here, [1b is Bohr's mag- 
neton. Evidently, the derived moments are larger than 
the localized spins of the individual Fe ions. They can be 
considered as the moments of small ferromagnetic clus- 
ters with perfectly parallel spins, in which the electrons 
can move freely, but with small misalignment between 



neighboring clusters. Since there are eight Fe 2+ ions per 
unit cell, their size ranges between about two and four 
unit cells. We also note that eq.(6) does not only de- 
scribe the magnetic field dependence of the MR effect 
but also nicely explains the decrease of the MR effect 
with increasing temperature. 

In summary, the MR effect of the Znj;Fc3_ x 04 films 
in the same way as the p X x(T) data can be well de- 
scribed within a model based on thermally activated hop- 
ping of itinerant charge carriers or small polarons be- 
tween the mixed- valent Fe^/Fe^ 4 " ions on the B sub- 
lattice. The activation energy contains a magnetic field 
dependent part, which depends on the misalignment 
of the neighboring spin moments. Within this model, 
MR^a; oc — 03(const. H/T) is expected in good agreement 
with the experimental data. In particular, the similar 
MK XX (H) curves for samples with strongly different p xx 
values and the decrease of the MR effect with increasing 
temperature are nicely reproduced. 



5. MAGNETIC PROPERTIES 

We have seen that Zn substitution and the creation 
of Fe vacancies in magnetite results in changes of the 
carrier density and hopping amplitudes of the itinerant 
charge carriers on the B sublattice. Following our dis- 
cussion in the introduction, this is expected to have sig- 
nificant influence on the magnetic properties. Therefore, 
we also systematically analyzed the magnetic properties 
of the Zn^Fea-zC^ films. The magnetic characteriza- 
tion was performed with the unpatterned films using 
SQUID magnetometry with magnetic fields hqH up to 
7T applied in the film plane. At room temperature, the 
M (H) loops show ferromagnetic behavior for all samples 
(Fig. 6). However, the measured saturation magnetiza- 
tion Ms as well as the remanent magnetization Mr and 
coercive field was found to strongly depend on both the 
Zn substitution level x and the growth atmosphere. 

We start our discussion with the stoichiometric FC3O4 
film grown in pure Ar atmosphere, serving as a refer- 
ence. Figure 6(a) shows that for this film the highest val- 
ues for the room temperature saturation magnetization 
(Ms = 3.16^B/f-u.) and remanence (Mr = 1.83 /is/f-u.) 
are obtained. We note that we do not obtain the theoret- 
ically expected value of Ms = 4//e/f-u. This is typical 
for thin film samples due to the presence of APBs 28,29 . 
However, the value measured for our film is among the 
highest reported in literature and corresponds well with 
the volume fraction of the APB phase of 20% estimated 
above. 

We next discuss the evolution of the saturation mag- 
netization of the Zn a; Fe3_ a ;04 films grown in pure Ar at- 
mosphere with increasing x. As shown in the inset of 
Fig. 6(a), both Ms and Mr were found to decrease with 
increasing x. Recalling eq.(l), the substitution of Fe^ + 
(3d 5 ,S = 5/2) by Zn 2 + (3d 10 , S = 0) results in a de- 
crease of the magnetization on the A sublattice. That 
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FIG. 6: (color online) Room temperature magnetization 
M versus magnetic field H applied in the film plane for 
Zn a; Fe3_ 2 ,04 thin films grown in (a) pure Ar atmosphere and 
(b) in an Ar/02 (99:1) mixture. The insets show the satura- 
tion magnetization Ms (full symbols) and the remanence Mr 
(open symbols) in units of ^/i.u. as a function of the Zn 
substitution level x. The lines are guides to the eye. 



not be directly compared. The reason is that the films in 
Ref. 12 have been grown at finite oxygen partial pressure 
of 1.0 x 10~ 6 mbar. In this case, the effect of Fc vacancies 
has to be taken into consideration as follows. 

The magnetic properties of the Zn a; Fe3_ 2: 04 films 
grown in an Ar/02 (99:1) mixture are shown in Fig. 6(b). 
There are pronounced differences to the films grown in 
pure Ar atmosphere. First, for x = both Ms and Mr 
are much lower than for the films grown in pure Ar. Sec- 
ond, Mg and Mr, is found to increase with increasing x 
for < x < 0.5 and then to decrease again for x = 0.9 
as shown in the inset of Fig. 6(b). The coercive field de- 
creases from about 55 mT to 8 mT on increasing x from 
to 0.9. The Ms(x) dependence can be straightforwardly 
explained by the presence of Fe vacancies in films grown 
in finite oxygen partial pressure and their removal by 
Zn substitution. For x = 0, the Fe vacancies shift the 
Fe^/Fe^ 1:1 balance towards Fe^ + , removing itinerant 
charge carriers and weakening the ferromagnetic double 
exchange at the expense of antifcrromagnctic superex- 
changc on the B sublatticc. This results in spin canting 
on the B sublattice, giving rise to a reduced saturation 
magnetization. On Zn substitution, part of the Fe vacan- 
cies are removed. This shifts the Fe^/Fe'g" ratio back 
towards a 1:1 balance. The related increase of the itin- 
erant charge carrier density and DE interaction on the 
B sublattice reduces the spin canting with increasing x. 
This explains the increase of the room temperature sat- 
uration magnetization in the range < x < 0.5 from 
M s = 1.05/i B /f.u. up to 3.02^ B /f-u. at x = 0.5. At 
10K, Ms increases from 1.26/UB/f-u. at x = up to 
3.92^ B /f-u. at x = 0.5. For even larger x, the Fe^ + 
moments on the A sublattice are strongly diluted, re- 
sulting in a strong weakening of the antiferromagnetic 



is, since the compensating magnetization of the A sub- 
lattice is removed, one would expect an increase of the 
total magnetization of the ferrimagnct. However, Zn sub- 
stitution also weakens Jab by diluting the A site mo- 
ments. Furthermore, it converts Fe 2 ^ (3d 6 ,S = 2) into 
Fe^ + on the B sublattice, resulting in a reduction of the 
itinerant charge carrier density. This weakens the ferro- 
magnetic double exchange, competing with antiferromag- 
netic superexchange interaction on the B sublattice. As 
discussed in the introduction, taken together this leads 
to an increase of the spin canting on the B sublattice 
with increasing x (cf. Fig. 1), explaining the observed 
reduction of the total magnetization of the ferrimagnetic 
Zn 2 ,Fc3_ a; 04 films. We note that recently Takaobushi 
et al. reported a monotonic increase of the saturation 
magnetization from about 0.5/ZB/f.u. to 3.2^iB/f-u. at 
10 K in the range < x < 0.9 and explained this ob- 
servation by the assumption that the A sublattice mag- 
netization is reduced by Zn substitution 12 . Comparison 
of the low temperature saturation magnetization of the 
x = films (3.31 /ie/f-u. when grown in Ar in our case, 
0.5 /ze/f.u. in Ref. 12 ) clearly shows that the samples can- 
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FIG. 7: (color online) Remanent magnetization Mr plotted 
versus temperature T for a Zno.5Fe2.5O4 film grown in an 
Ar/02 (99:1) mixture. The data were taken after field cooling 
at 7 T on increasing the temperature at zero applied magnetic 
field. The inset schematically shows the hump in the M(T) 
curve as the result of two opposite sublattice magnetizations 
with different T dependence. 
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superexchange Jab between the A and B sublattices. 
Furthermore, Zn substitution may overcompensate the 
Fe vacancies again shifting the Fcg + /Fc^ + balance to- 
wards Fe^ + . Taken together, this explains the decrease 
of Ms above the critical value of x f» 0.5, which of course 
depends on the oxygen partial pressure during growth 
and the corresponding amount of Fe vacancies. The ob- 
served increase of Ms with increasing x for Zn a ,Fe3_ a; 04 
films grown in an Ar/02 mixture is consistent with the 
result of Takaobushi et al. 12 , who reported an increase 
of Mg(lOK) from 0.5/is/f-u. to 3.2/is/f-u. on increas- 
ing x from to even 0.9. However, our systematic study 
strongly suggests that this increase is most likely caused 
by the removal of Fe vacancies rather than due to the 
removal of Fe^ + moments on the A sites as argued in 
Ref. 12 . 

Figure 7 shows the temperature dependence of the re- 
manent magnetization Mr for a Zno.5Fe2.5O4 film grown 
in an Ar/02 mixture. In agreement with the p xx (T) data 
of Fig. 4, there is no indication for the Verwey transition 
in Mr(T) around 120 K. This is the case for all samples, 
except for the FC3O4 film grown in pure Ar. This obser- 
vation is not only consistent with our transport data but 
also with magnetization data from other groups, indicat- 
ing a smearing or suppression of the Verwey transition 
by an amount of Zn substitution as small 0.01 68 . 
As shown in Fig. 7, the measured Mr (T) curve first de- 
creases with increasing T, goes through a minimum at 
T rj 150 K, then increases up to T w 275 K, before it 
decreases again towards room temperature. Such a be- 
havior is often observed in fcrrimagnets and can be ex- 
plained on the basis of Neel's two sublatticc model 69 . 
Within this model the origin of the hump at 275 K is 
a different T dependence of the two sublattice magne- 
tizations Ma{T) and Mb(T) due to different effective 
molecular fields on the A and B sites. In our case the ef- 
fective molecular field on the A sublattice is expected to 
be smaller due to the partial substitution of Fe^ + by non- 
magnetic Zn 2+ . As a result, the combined magnetization 
M(T) = M A (T)+M B {T) should show an upward hump 
in agreement with the experimental result. Note that 
the observation of a hump, which is directed upwards, 
provides direct experimental evidence for A site substi- 
tution of Zn 2+ in the inverse spinel structure 69 . For the 
x = 0.33 sample, no hump could be observed, most likely 
due to a too small imbalance in the effective molecular 
fields at the smaller substitution level. 

In conclusion, the magnetic properties of the 
Zn :c Fe3_ x 04 films grown under different oxygen par- 
tial pressure strongly depend on whether the films were 
grown in pure Ar atmosphere or at finite oxygen par- 
tial pressure. Fully consistent with the transport data 
these differences can be explained by the following ef- 
fects: (i) For films grown in pure Ar, Zn substitution 
on the A site removes Fe^ + moments from the A sublat- 
tice and due to charge neutrality itinerant electrons from 
the B sublattice. The resulting weakening of both the 
antifcrromagnctic A-O-B SE and the ferromagnetic DE 



1 1 1 1 1 1 1 1 1 1 II 

: Zn Fe, x = 

x 3-x 4 _ 

f 7=300K / 05 

/ 


1 1 1 1 1 1 
x = 

W 0.0 - 

/ - 

0.1 


~ • n **** 

■ 0.5 


grown in - 


'_ / 

/ 

/ 


nurp Ar 

uui ^. ni 


/ 

- ' grown in 




- / Ar/O. mixture 

/ 2 




'- *0.0 

1 1 1 II 


, , , , I , " 



10 2 10 3 10 4 

a (Q" 1 m 1 ) 

XX x ' 



FIG. 8: (color online) Saturation magnetization Ms plot- 
ted versus the longitudinal conductivity a xx for epitaxial 
Zna;Fe3_ x 04 thin films with different Zn substitution x < 0.5 
grown in pure Ar atmosphere (squares) or an Ar/C>2 (99:1) 
mixture (circles). 



on the B sublattice results in an increasing spin canting 
on the B sublattice with increasing x. This causes a re- 
duction of the saturation magnetization with increasing 
x. (ii) For films grown in oxygen partial pressure, the 
magnetic properties are strongly influenced by the pres- 
ence of Fe vacancies. They reduce the density of itinerant 
electrons on the B sublattice thereby weakening the fer- 
romagnetic DE at the expense of antiferromagnetic SE. 
This causes spin canting on the B sublattice and, in turn, 
a reduced saturation magnetization. The effect of addi- 
tional Zn substitution is the partial removal of the Fe 
vacancies. This effectively increases the density of itiner- 
ant electrons on the B sublattice thereby strengthening 
the ferromagnetic DE. The result is an increase of the 
saturation magnetization with increasing x. Again, the 
observed dependence of the saturation magnetization on 
the Zn substitution levels x cannot be simply explained 
by different APB densities in the different thin film sam- 
ples. It is very unlikely that the APB density increases 
with increasing x for samples of the same thickness when 
grown in Ar atmosphere while for those prepared in the 
Ar/02 mixture it first decreases for x < 0.5 and then 
increases again for x > 0.5. 

We finally link the magnetic properties to the transport 
data by plotting the saturation magnetization Ms versus 
the electrical conductivity a xx ~ l/p xx 44 for films with 
various x and grown in Ar or Ar/02 (Fig. 8). Evidently, 
there is a clear correlation between Mg and a xx . Samples 
with high Ms have large a xx and vice versa. Again, this 
correlation can be straightforwardly explained within the 
scenario presented above. Charge transport is dominated 
by the hopping of the itinerant ti g electrons between the 
mixed- valcnt Fc^/Fc^ ions on the B sublatticc. Due 
to the strong on-site Hund's rule coupling the spin of the 
hopping spin-down electron is aligned anti-parallel to the 
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local moment of the spin- up electrons (cf. Fig. 1). There- 
fore, the hopping amplitude is significantly suppressed if 
the moments of neighboring B sites are not parallel. That 
is, any spin canting on the B sublattice reducing the sat- 
uration magnetization also results in a reduction of the 
electrical conductivity. 

6. CONCLUSION 

We have grown epitaxial Zn x Fc3_ x 04 thin films (0 < 
x < 0.9) with excellent structural properties both in pure 
Ar atmosphere and an Ar/02 mixture using laser molecu- 
lar beam epitaxy. We show that the electrical conductiv- 
ity and the saturation magnetization can be tuned over a 
wide range (1.0 ... 3.2 ^ B /f-u. and 10 2 . . . 10 4 n^m -1 at 
room temperature) by Zn substitution and/or finite oxy- 
gen partial pressure during growth. This demonstrates 
that Zn x Fc3_ x 04 can be used as a material with tailored 
electrical and magnetic properties. Our comprehensive 
study shows that electrical conduction, magnetotrans- 
port, and the magnetic properties of Zn 2 ,Fe3_ a; 04 are 
dominated by the density and thermally activated hop- 
ping of the itinerant t2 g electrons between the mixcd- 
valent Fe| f /Fe| + ions, providing electrical conduction 
and mediating a ferromagnetic DE interaction on the B 
sublattice. This scenario is confirmed by the observation 
that the saturation magnetization scales with the longi- 
tudinal conductivity and the field dependence of the MR 
effect. Our combined systematic analysis of the mag- 
netotransport and magnetic properties also provides a 
clear picture of the effect of Zn substitution and growth 
in finite oxygen partial pressure. First, Zn substitution 



removes both Fe^ moments from the A sublattice and 
itinerant charge carriers from the B sublattice. Second, 
growth in finite oxygen partial pressure generates Fe va- 
cancies, also removing itinerant charge carriers from the 
B sublattice. Hence, in both cases a reduction of the 
charge carrier density and a weakening of the ferromag- 
netic DE on the B sublattice is obtained. This results in 
an increase of resistivity and a reduction of the saturation 
magnetization due to spin canting on the B sublattice. 
Third, applying both Zn substitution and growth in oxy- 
gen at the same time does not result in an additive effect. 
In contrast, a compensation effect is observed, since Zn 
substitution removes part of the Fe vacancies. To sum 
up, one can say that the high Curie temperature together 
with the electrical and magnetic properties that can be 
tuned over a wide range make Zn x Fc3_ T 04 a promising 
functional material for spintronic devices. The combina- 
tion of Zn 2; Fc3_ : j;04 with semiconductor materials such 
as ZnO in multi-functional heterostructures 27 will be par- 
ticularly appealing. 

Acknowledgments 



We thank Andreas Erb for the preparation of poly- 
crystalline target materials for the laser-MBE process 
and Thomas Brcnningcr for continuous technical sup- 
port. Financial support by the German Science Founda- 
tion within the priority programs 1157 and 1285 (project 
Nos. GR 1132/13 & 14) and the German Excellence Ini- 
tiative via the Nanosystems Initiative Munich (NIM) is 
gratefully acknowledged. D.V. and M.S.R.R. thank the 
DA AD for financial support. 



* Electronic address: Matthias.Opel@wmi.badw.de 
t Electronic address: Rudolf.Gross@wmi.badw.de 

1 S.A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. 
Daughton, S. von Molnar, M. L. Roukes, A. Y. Chtchelka- 
nova, and D. M. Treger, Science 294, 1488 (2001). 

2 I. Zutic, J. Fabian, S. Das Sarma, Rev. Mod. Phys. 76, 
323 (2004). 

3 J. S. Moodera, L. R. Kinder, T. M. Wong, R. Meservey, 
Phys. Rev. Lett. 74, 3273 (1995). 

4 H. Ohno, D. Chiba, F. Matsukura, T. Omiya, E. Abe, T. 
Dietl, Y. Ohno, and K. Ohtani, Nature 408, 944 (2000). 

5 D. Chiba, M. Yamanouchi, F. Matsukura, and H. Ohno, 
Science 301, 943 (2003). 

6 S.T.B. Goennenwein, M. Althammer, C. Bihler, A. Brandl- 
maier, S. Geprags, M. Opel, R. Gross, W. Schoch, W. Lim- 
mer, H. Huebl, and M. S. Brandt, phys. stat. sol. (RRL) 
2, 96 (2008). 

7 C. Bihler, M. Althammer, A. Brandlmaier, S. Geprags, 
M. Weiler, M. Opel, W. Schoch, W. Limmer, R. Gross, 
M. S. Brandt, and S.T.B. Goennenwein, Phys. Rev. B 78, 
045203 (2008). 

8 H. Tanaka, J. Zhang, and T. Kawai, Phys. Rev. Lett. 88, 
027204 (2001). 



9 Z. Zhang and S. Satpathy, Phys. Rev. B 44, 13319 (1991). 

10 Y. S. Dedkov, U. Riidiger, and G. Giintherodt, Phys. Rev. 
B 65, 064417 (2002). 

11 M. Fonin, Y. S. Dedkov, R. Pentcheva, U. Riidiger and G. 
Giintherodt, J. Phys.: Condens. Matter 19, 315217 (2007); 
see also J. Phys.: Condens. Matter 20, 142201 (2008). 

12 J. Takaobushi, H. Tanaka, T. Kawai, S. Ueda, J.-J. Kim, 
M. Kobata, E. Ikenaga, M. Yabashi, K. Kobayashi, Y. 
Nishino, D. Miwa, K. Tamasaku, and T. Ishikawa, Appl. 
Phys. Lett. 89, 242507 (2006). 

13 J. Takaobushi, M. Ishikawa, S. Ueda, E. Ikenaga, J.-J. 
Kim, M. Kobata, Y. Takeda, Y. Saitoh, M. Yabashi, Y. 
Nishino, D. Miwa, K. Tamasaku, T. Ishikawa, I. Satoh, H. 
Tanaka, K. Kobayashi, and T. Kawai, Phys. Rev. B 76, 
205108 (2007). 

14 A. Brandlmaier, S. Geprags, M. Weiler, A. Boger, M. 
Opel, H. Huebl, C. Bihler, M.S. Brandt, B. Botters, D. 
Gmndler, R. Gross, and S.T.B. Goennenwein, Phys. Rev. 
B 77, 104445 (2008). 

15 L. Neel, Ann. Phys. 3, 137 (1948). 

16 Y. Yafet and C. Kittel, Phys. Rev. 87, 290 (1952). 

17 J. Loos and P. Novak, Phys. Rev. B 66, 132403 (2002) 

18 R. J. McQueeney, M. Yethiraj, S. Chang, W. Montfrooij, 



13 



T. G. Perring, J. M. Honig, and P. Metcalf, Phys. Rev. 
Lett. 99, 246401 (2007). 

19 A. Rosencwaig, Phys. Rev. 181, 946 (1969). 

20 P. Wang, Z. Kakol, M. Wittenauer, and J. M. Honig, Phys. 
Rev. B 42, 4553 (1990). 

21 V. G. Harris, N. C. Koon, C. M. Williams, Q. Zhang, M. 
Abe, and J. P. Kirkland, Appl. Phys. Lett. 68, 2082 (1996). 

22 Y. Li, Q. Li, M. Wen, Y. Zhang, Y. Zhai, Z. Xie, F. Xu, 
and S. Wei, J. Elect. Spectrosc. 160, 1 (2007). 

23 J.P. Shepherd, R. Aragon, J.W. Koenitzer, and J.M. 
Honig, Phys. Rev. B 32, 1818 (1985). 

24 R. Aragon and R. H. McCallister, Phys. Chem. Minerals 
8 112 (1982). 

25 R. Aragon, D.J. Buttrey, J.P. Shepherd, and J.M. Honig, 
Phys. Rev. B 31, 430 (1985). 

26 K. A. Shaw, E. Lochner, and D. M. Lind, J. Appl. Phys. 
87, 1727 (2000). 

27 A. Nielsen, A. Brandlmaier, M. Althammer, W. Kaiser, M. 
Opel, J. Simon, W. Mader, S.T.B. Goennenwein, and R. 
Gross, Appl. Phys. Lett. 93, 162510 (2008). 

28 D. T. Margulies, F. T. Parker, F. E. Spada, R. S. Goldman, 
J. Li, R. Sinclair, and A. E. Berkowitz, Phys. Rev. B 53, 
9175 (1996); see also D. T. Margulies, F. T. Parker, M. L. 
Rudee, F. E. Spada, J. N. Chapman, P. R. Aitchison, and 
A. E. Berkowitz, Phys. Rev. Lett. 79, 5162 (1997). 

29 T. Hibma, F. C. Voogt, L. Niesen, P. A. A. van der Heijden, 
W. J. M. de Jonge, J. J. T. M. Donkers, and P. J. van der 
Zaag, J. Appl. Phys. 85, 5291 (1999). 

30 R. Gross, J. Klein, B. Wiedenhorst, C. Hofener, U. Schoop, 
J. B. Philipp, M. Schonecke, F. Herbstritt, L. Alff, Yafeng 
Lu, A. Marx, S. Schymon, S. Thienhaus, and W. Mader, 
Proc. SPIE 4058, 278 (2000). 

31 M.E. Fleet, Acta Cryst. B 37, 917 (1981). 

32 J. Klein, C. Hofener, L. Alff, and R. Gross, Supercond. Sci. 
Technol. 12, 1023 (1999). 

33 D. Reisinger, B. Blass, J. Klein, J. B. Philipp, M. Scho- 
necke, A. Erb, L. Alff, and R. Gross, Appl. Phys. A 77, 
619 (2003). 

34 D. Reisinger, M. Schonecke, T. Brenninger, M. Opel, A. 
Erb, L. Alff, R. Gross, J. Appl. Phys. 94 1857 (2003). 

35 D. Reisinger, P. Majewski, M. Opel, L. Alff, R. Gross, 
Appl. Phys. Lett. 85, 4980 (2004). 

36 F. Schedin, E. W. Hill, G. van der Laan, and G. Thornton, 
J. Appl. Phys. 96, 1165 (2004). 

37 A. Marcu, T. Yanagida, K. Nagashima, H. Tanaka, and T. 
Kawai, J. Appl. Phys. 102, 023713 (2007). 

38 I. Barin, Thermodynamical Data of Pure Substance, British 
Library Cataloguing in Publication Data, VCH Verlagsge- 
sellschaft, Weinheim (1989). 

39 I.G. Austin and N.F. Mott, Adv. Phys. 18, 41 (1969). 

40 N.F. Mott and Z. Zinamon, Rep. Prog. Phys. 33, 881 
(1970). 

41 D. Emin and T. Holstein, Ann. Phys. 53, 439 (1969). 

42 W. Haubenreisser, phys. stat. sol. 1, 619 (1961). 

43 L. Alff, G. Fischer, R. Gross, F. Kober, K.D. Husemann, 
A. Beck, T. Nissel, C. Burckhardt, and F. Schmidl, Physica 



C 200, 277 (1992). 

44 D. Venkateshvaran, W. Kaiser, A. Boger, M. Althammer, 
M.S. Ramachandra Rao, S.T.B. Goennenwein, M. Opel, 
R. Gross, Phys. Rev. B 78, 092405 (2008). 

45 A.J. Bosman and H.J. van Daal, Adv. Phys. 19, 1 (1970). 

46 M.I. Klinger and A. A. Samokhvalov, phys. stat. sol. (b) 
79, 9 (1977). 

47 N. Popandian and A. Narayanasamy, J. Appl. Phys. 93, 
2770 (2002). 

48 Z. Simsa and V.A.M. Brabers, IEEE Trans. Magn. 24, 
1910 (1988). 

49 G. J. Snyder, R. Hiskes, S. DiCarolis, M. R. Beasley, and 
T. H. Geballe, Phys. Rev. B 53, 14434 (1996). 

50 Yafeng Lu, J. Klein, F. Herbstritt, J. B. Philipp, A. Marx, 
and R. Gross, Phys. Rev. B 73, 184406 (2006). 

51 Yafeng Lu, J. Klein, F. Herbstritt, J.B. Philipp, A. Marx, 
L. Alff, and R. Gross, phys. stat. sol. (b) 242, 1545 (2005). 

52 Yafeng Lu, J. Klein, C. Hofener, B. Wiedenhorst, J. B. 
Philipp, F. Herbstritt, A. Marx, L. Alff, and R. Gross, 
Phys. Rev B 62, 15806 (2000). 

53 E. J. W. Verwey, Nature (London) 144, 327 (1939). 

54 N. F. Mott, Adv. Phys. 16, 49 (1967). 

55 E. J. Verwey, P.W. Haayman, and F. C. Romeijn, J. Chem. 
Phys. 15, 181 (1947). 

56 J. Garcia, G. Subfas, M. G. Proietti, J. Blasco, H. 
Renevier, J. L. Hodeau, and Y. Joly, Phys. Rev. B 63, 
054110 (2001); see also J. Garcia, G. Subfas, M. G. Proi- 
etti, H. Renevier, Y. Joly, J. L. Hodeau, J. Blasco, M. C. 
Sanchez, and J. F.Berar, Phys. Rev. Lett. 85, 578 (2000). 

57 F. Walz, J. Phys.: Condens. Matter 14, R285-R340 (2002). 

58 J. P. Wright, J. P. Attfield, and P. G. Radaelli, Phys. Rev. 
B 66, 214422 (2002). 

59 I. Leonov, A. N. Yaresko, V. N. Antonov, M. A. Korotin, 
and V. I. Anisimov, Phys. Rev. Lett. 93, 146404 (2004). 

60 S. Todo, K. Siratori, and S. Kimura, J. Phys. Soc. Japan 
64, 2118 (1995). 

61 W. Eerenstein, T. T. M. Palstra, and T. Hibma, Phys. 
Rev. B 66, 201101 (2002). 

62 S. Tiwari, D. M. Phase, and R. J. Choudharya, Appl. Phys. 
Lett. 93, 234108 (2008). 

63 M. L. Parames, J. Mariano, Z. Viskadourakis, N. Popovici, 
M. S. Rogalski, J. Giapintzakis, and O. Conde, Appl. Surf. 
Sci. 252, 4610 (2006). 

64 G. A. Ferguson jr and M. Hass, Phys. Rev. 112, 1130 
(1958). 

65 M. Viret and L. Ranno, J.M.D. Coey, Phys. Rev. B 55, 
8067 (1997). 

66 M. Viret, L. Ranno, and J.M.D. Coey, J. Appl. Phys. 81, 
(1997). 

67 P. Wagner, I. Gordon, L. Trappeniers, J. Vanacken, F. 
Herlach, V. V. Moshchalkov, and Y. Bruynseraede , Phys. 
Rev. Lett. 81, 3980 (1998). 

68 J. M. de Teresa, private communication (2008). 

69 A. H. Morrish, The physical principles of magnetism, IEEE 
Press, New York (2001). 



